The twilight zone in the parametric evolution of eigenstates: beyond perturbation theory and 

semiclassics 
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Considering a quantized chaotic system, we analyze the evolution of its eigenstates as a result of varying a 
control parameter. As the induced perturbation becomes larger, there is a crossover from a perturbative to a 
non-perturbative regime, which is reflected in the structural changes of the local density of states. For the first 
time the full scenario is explored for a physical system: an Aharonov-Bohm cylindrical billiard. As we vary 
the magnetic flux, we discover an intermediate twilight regime where perturbative and semiclassical features 
co-exist. This is in contrast with the simple crossover from a Lorentzian to a semicircle line-shape which is 
found in random-matrix models. 
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The analysis of the evolution of eigenvalues and of the 
structural changes that the corresponding eigenstates of a 
chaotic system exhibit as one varies a parameter cj) of the 
Hamiltonian 7i((/)) has sparked a great deal of research ac- 
tivity for many years. Physically the change of (j> may repre- 
sent the effect of some externally controlled field (like electric 
field, magnetic flux, gate voltage) or a change of an effective- 
interaction (as in molecular dynamics). Thus, these studies 
are relevant for diverse areas of physics ranging from nuclear 
|Ql,i2i] and atomic physics L3,,4J to quantum chaos liSn^iTuSJ 
and mesoscopics I9tll0ll. 

Up to now the majority of this research activity was focused 
on the study of eigenvalues, where a good understanding has 
been achieved, while much less is known about eigenstates. 
The pioneering work in this field has been done by Wigner 
II2I], who studied the parametric evolution of eigenstates of a 
simplified Random Matrix Theory (RMT) model of the type 
H = E + (pB. The elements of the diagonal matrix E are the 
ordered energies {En}, with mean level spacing A, while B 
is a banded random matrix. Wigner found that as the param- 
eter (j> increases the eigenstates undergoes a transition from a 
perturbative Lorentzian-type line shape to a non-perturbative 
semicircle line-shape. 

For many years the study of parametric evolution for canon- 
ically quantized systems was restricted to the exploration of 
the crossover from integrability to chaos [7, Only later 
0] it has been realized that a theory is lacking for systems 
that are chaotic to begin with. Inspired by Wigner theory, the 
natural prediction was that the local density of states (LDOS) 
should exhibit a crossover from a regime where a perturbative 
treatment is applicable, to a regime where semiclassical ap- 
proximation is valid. However, despite a considerable amount 
of numerical efforts 0], there was no clear-cut demonstration 
of this crossover Neither a theory has been developed de- 
scribing how the transition from the perturbative to the non- 
perturbative regime takes place. 

It is the purpose of this Letter to present, for the first time, a 
complete scenario of parametric evolution, in case of a phys- 
ical system that exhibits hard chaos. We explore the validity 
of perturbation theory and semiclassics, and we discover the 
appearance of an intermediate regime ("twilight zone") where 



both perturbative and semiclassical features co-exist. With- 
out loss of generality we consider as an example a billiard 
system whose classical dynamics is characterized by a corre- 
lation time Ta, which is simply the ballistic time. Associated 
with r,i is the energy scale h/xa. Next we look on a similar 
billiard, but with a rough boundary. This roughness is charac- 
terized by a length scale which is £ times smaller, hence we 
can associate with it an energy scale JE'nu = (^/''"d) ^- The 
roughness does not affect the chaoticity: the coirelation time 
Ta as well as the whole power spectrum are barely affected. 
Consequently we explain that (5i?Nu is not reflected in the RMT 
modeling of the Hamiltonian. Still in the LDOS analysis we 
find that non-universal (system specific) features appear. The 
appearance of such features is a generic phenomenon in quan- 
tum chaos studies. It introduces a new ingredient into the the- 
ory of parametric evolution which goes beyond RMT. 

The model that we will use in our analysis is a particle 
confined to an Aharonov-Bohm (AB) cylindrical billiard (see 
Fig.l) where one can control the magnetic flux $. The cylin- 
drical billiard is constructed by wrapping a 2D billiard with 
hard wall boundaries. The lower boundary at y = is flat, 
while the upper boundary y — Ly -\- W^{x) is deformed. The 
deformation is described by S,{x) = X]ri=i cos(na;) where 
An are random numbers in the range [-1,1]. The illustration in 
Fig.l assumes a smooth boundary (i — 1). The Hamiltonian 
of a particle in the cylindrical AB billiard is 
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y = L + W^(x) 




FIG. 1: Left: 2D billiard with 
Aharonov-Bohm cylinder. 



1. Right: Corresponding 
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supplemented by periodic boundary conditions in the hor- 
izontal direction, and hard wall boundary conditions along the 
lower and upper boundaries, and py are the momenta. 
Later we shall use the notation = e^/h. We consider the 
chaotic 7i((/i = 0) as the unperturbed Hamiltonian. 

After conformal transformation |7] the billiard is mapped 
into a rectangular, with a mass tensor which is space depen- 
dent. Then it is possible to compute the matrix representation 
of the Hamiltonian in the plane wave basis ji^/i) of the rectan- 
gular. The result is: 



t(0,2) , 2 .(2,2)4 TT^/^'- 



2^(j£+J^ (2,2) . ly+l^' - <j)/n (1,1) 
e , „ .-.2 '^i^'i' *^ ,,2 _ ,,/2 "^i^'i^ 



(2) 



where 



j(i,k) _ 



The classical dimensionless parameters of the model are the 
aspect ratio a = Ly/Lx, the tilt relative amplitude e = 
W/Ly, and the roughness parameter I. Upon quantization 
we have h that together with m and E determines the De- 
Broglie wavelength of the particle, and hence leads to an ad- 
ditional dimensionless parameter [LxLy/ {2'Kh'^)]mE. 
For 2D billiards the mean level spacing A is constant, and 
hence ue — E/ A cx can be interpreted as either the 
scaled energy or as the level index. Optionally we define a 
semiclassical parameter /ijcaicd = l/^/n^. 

In the numerical study we have taken e=0.06 and a—1, 
for which the classical dynamics is completely chaotic (for 
any (j)). We consider either £=1 for smooth boundary, or 
£=100 for rough boundary. The eigenstates of the 

Hamiltonian 7i(0) were found numerically for various val- 
ues of the flux (0.0006 < (j> < 60). We were interested in 
the states within an energy window SE^Ab that contains 
i5n£;~200 levels around the energy i?«400. Note that the 
size of the energy window is classically small (5E ^ E), but 
quantum mechanically large (SE ^ A). 

The object of our interest are the overlaps of the eigen- 
states \n{(/))) with a given eigenstate |m(0)) of the unper- 
turbed Hamiltonian: 



P(n|m) = |(n(</))|m(0))p = 



dxdydpxdpy_^^^)^^^) 



The overlaps P{n\m) can be regarded as a distribution with 
respect to n. Up to some trivial scaling it is essentially the 
local density of states (LDOS). The associated dispersion is 
defined as 6E = P[n\m){En — EmYY^'^ practice we 
plot P{n\m) as a function of r = n — m or as a function 
of {En— Em), and average over the reference state m. The 



second equality in Q is useful for the semiclassical analy- 
sis. It involves the Wigner functions p^^^\x, y,Px,Py) which 
are associated with the eigenstates |n(</))). The semiclassical 
approximation is based on the microcanonical approximation 
p^"^ c>c 6{En—Ti.{x,y,Px,Py))- With this approximation the 
integral can be calculated analytically leading to 



Pc,(n\m) 



(4) 



7rJ2{SE,,f - [{En~Em) - 5Ell{2E„ 



where SE^ — {hvE/Lx)4> with ve — (2i?/m)^/^. It is im- 
plicit in that Pa{n\m)—Q outside of the allowed range, 
which is where the expression under the square root is neg- 
ative: For large \En~E„i\ there is no intersection of the cor- 
responding energy surfaces, and hence no classical overlap. 

A few words are in order regarding quantum to classical 
correspondence (QCC). Whenever P{n\m) sa P^f{n\m) we 
call it "detailed QCC", while 5E w 5E^^ is referred to as "re- 
stricted QCC" 1 6]. It is remarkable that (the robust) restricted 
QCC holds even if (the fragile) detailed QCC fails completely. 
We have verified 1 13] that also in the present system 5E is nu- 
merically indistinguishable from SE^- 

A fixed assumption of this work is that is classically 
small. But quantum mechanically it can be either 'small' or 
'large'. Quantum mechanically small (p means that pertur- 
bation theory do provide a valid approximation for P{n\m). 
What is the border between the perturbative regime and the 
non-perturbative regime, we discuss later First we would 
like to show that the prediction which is based on perturba- 
tion theory, to be denoted as Pp„(n|m), is very different from 
the semiclassical approximation. 

In order to write the expression for Pp„(n|m) we have first 
to clarify how to apply perturbation theory in the context of 
the present model. To this end, we write the perturbed Hamil- 
tonian 7i(0) in the basis of 7i((/) — 0). Since we assume that 
the perturbation is classically small, it follows that we can lin- 
earize the Hamiltonian with respect to Consequently the 
perturbed Hamiltonian is written as 7i = E+(f>B, where E = 
diag{i?„} is a diagonal matrix, while B = { — (^/e)X,„„}. 
The current operator is conventionally defined as 

I=-dn/d^^{e/imLx))px 

Its matrix elements can be found using a semiclassical recipe 
IIU, namely |I„,„|2 ^ {A/{2TTh))CiiEn-E„,)/h), where 
C{uj) is the Fourier transform of the current-current corre- 
lation function C(r). Conventional condensed matter cal- 
culations are done for disordered rings where one assumes 
C(r) to be exponential, with time constant r^, which is essen- 
tially the balUstic time. Hence C{u;) oc l/(w^ + (1/rd)^) 
is a Lorentzian. This Lorentzian approximation works well 
also for the chaotic ring that we consider In fact we can 
do better by exploiting a relation between T{t) and the force 
T{t) = ~px, leading to Ciuj) = {e/(mLx)fC:F[^)/uj^. 
The force !F{t) is a train of spikes corresponding to collisions 
with the boundaries. Assuming that the collisions are uncor- 
rected on short times we have Cjr{uo) 



i8/3n)m\%/Ly, 
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FIG. 2: (a) The classical power spectram C{uj) plotted 
(in grey) together with the quantum mechanical band-profile 
(27reVAft)|B„„|2 fori = 1, and /iscaied 0.018. (b) The LDOS 
kernel P{n\m) in the perturbative regime for a billiard with £ = 100 
and perturbation (p — 2.7. (c) Same as (b) but zoomed normal scale. 
The width of the non-perturbative component is F/A — 36. Note 
that in this regime the variance AE/ A ~ 58 is still dominated by 
the (perturbative) tails. For comparison we display the calculated 
Pel, and Prmt. 



> (l/Td). This is known as the "white noise" approx- 



for Lu ^ 

imation I12ll . We have checked the vaUdity of this approxi- 
mation in the present context by a direct numerical evaluation 
of C{uj), and also verified the validity of the above recipe by 
direct evaluation of the matrix elements of B via Eq.lEli, see 
Fig. 2(a). The classical C{uj) was numerically evaluated by 
Fourier analysis of the fluctuating current T{t) for a very long 
ergodic trajectory that covers densely the whole energy sur- 
face H{0) = E. 

Perturbation theory to infinite order with the Hamiltonian 
H ~ E + (f)B leads to a Lorentzian-type approximation 
for the LDOS |2] (see also Section 18 of |6']c). It is an ap- 
proximation because all the higher orders are treated within 
a Markovian-like approach (by iterating the first order result) 
and convergence of the expansion is pre-assumed, leading to 
Pp„(n|m) = (l)^\Bn„,\^ /[T^ + {En-E„,)% In practice the pa- 
rameter r{(j)) can be determined (for a given 0) by imposing 
the requirement of having -Pprt(r) normalized to unity. Substi- 
tuting the expression for the matrix elements we get 



P^„{n\m) 



Sh\hvEf/{i^^LlLl) 



(5) 



By comparing the exact P{r) to the approximation Eq.(|5} we 
can determine the regime (f) < (p^^ for which the approxima- 
tion P{r) Ri Ppnir) makes sense. The practical procedure to 
determine 0p„ is to plot SE^^ and to see where it departs from 
SEc,. The latter is a linear function of (p while the former be- 
comes sublinear for large enough (j>, (and even would exhibit 
saturation if we had a finite bandwidth). In case of Eq.(|5|i this 
reasoning leads to a crossover when SEa{(f>) ^ fi/ra. Hence 
we g et that the border of the perturbative regime (see footnote 

) is (ppn = La;/ (vETa) ~ 1. 

What happens to P{r) in practice? If we take the Wigner 
RMT model as an inspiration, we expect to have at (p ^ (p^,, 
a simple crossover from a Pp,, line-shape to a P, line-shape. 
The latter is regarded as the semiclassical analogue of the (ar- 
tificial) semicircle line shape. Indeed for the smooth billiard 
(£ = 1) we have verified that this naive expectation is realized 
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FIG. 3: (a) The LDOS kernel P(n]m) for (?!) = 31.4, where ^= 100. 
(b) The same parameters but £ = 1. In panel (a) we observe co- 
existence of perturbative and SC structures while in panel (b) we 
witness detailed QCC. 



fv^. But for the rough billiard {I = 100) we witness a more 
complicated scenario. In Fig. 2(b,c) we show the LDOS for 
<p < 0p„, where it (still) agrees quite well with Pp„. In Fig. 3 we 
show the LDOS for (p > 0p,i, where we would naively expect 
agreement with P,,. Rather we witness a three peak structure, 
where the r ~ peak is of perturbative nature, while the other 
are the fingerprint of semiclassics. For sake of comparison we 
show the corresponding results for a smooth billiard (i = 1) 
and otherwise the same parameters. There we have detailed 
QCC as is naively expected. The co-existence of perturba- 
tive and semiclassical features persists within an intermediate 
regime of (p values, to which we refer as the "twilight zone". 

Before we adopt a phase space picture in order to explain 
the above observations, we would like to verify that indeed 
random matrix modeling does not lead to a similar effect: Af- 
ter all the standard Wigner model, that gives rise to a simple 
crossover from a Lorentzian to a semicircle line shape, as- 
sumes a simple banded matrix, which is not the case in our 
model. As argued above the matrix elements of B decay as 
l/\n — mp from the diagonal. This implies that Ppri(?') is in 
fact not a Lorentzian, and also may imply that the crossover 
to the non-perturbative regime is more complicated. In order 
to resolve this subtlety we have taken a randomized version of 
the Hamiltonian H = E + (pB. Namely, we have random- 
ized the signs of the off-diagonal elements of the B matrix. 
Thus we get an RMT model with the same band profile as 
in the physical model. This means that Pp„ is the same for 
both models (the physical and the randomized), but still they 
can differ in the non-perturbative regime. Indeed, looking at 
the LDOS of the randomized model we observe that the semi- 
classical features are absent: Prmt(7') unlike P(r) exhibits a 
simple crossover from perturbative to non-perturbative line- 
shape. 

In what follows we would like to argue that the structure 
of P(r), both perturbative and non-perturbative components, 
can be explained using a phase space picture. [For phrasing 
purpose we find the "Wigner function language" most conve- 
nient, still the reader should notice that we do not need or use 
this representation in practice]. We recall the P{n\m) is deter- 
mined by the overlap of two Wigner functions. In the present 
context the Wigner functions p'") are supported by shifted cir- 
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FIG. 4: The probability distribution \{v,^\n)\^ for (a) The n = 
2423 eigenstate of the smooth {t = 1) billiard; (b) The n = 1000 
eigenstate of the rough (I = 100) billiard. Note that this is essentially 
the (j>x,Py) momentum distribution. The state in panel (a), unlike 
the state in panel (b), is a typical semiclassical state. Namely it is 
well concentrated on the energy shell. 



cles {px — {(j) / {2tt)Y +Py — 2mEn- We are looking for their 
overlap with a reference Wigner function which is supported 
by the circle p'^ + Py = 2m Em. The question is whether 
the overlaps of the Wigner functions p*^"^ and p'™^ can be 
approximated by a classical calculation, and under what cir- 
cumstances we need perturbation theory. 

Generically the Wigner function has a transverse Airy-type 
structure. If the "thickness" of the Wigner function is much 
smaller compared with the separation \En—Em\ of the en- 
ergy surfaces then we can trust the semiclassical approxima- 
tion. This will always be the case if h is small enough, or 
equivalently if we can make <j) large enough. In such case the 
dominant contribution comes from the intersection of the en- 
ergy surfaces, which is the phase space analogue of stationary 
phase approximation. The other extreme is the case where the 
"thickness" of Wigner function is larger compared with the 
separation of the energy surfaces (namely 5Ea{<t>) < h/T^). 
Then the contribution to the overlap comes "collectively" 
from all the regions of the Wigner (quasi) distribution, not 
just from the intersections. In such case we expect perturba- 
tion theory to work. 

The above reasoning assumes that the wavefunction is con- 
centrated in an ergodic-like fashion in the vicinity of the en- 
ergy surface. This is known as "Berry conjecture" fTTj. In 
case of billiards it implies that the wavefunction looks like a 
random superposition of plane waves with \p\ = (2mi?)^/^. 
We find (see Fig. 4) that this does not hold in case of a rough 
billiard (unless h were extremely small, so as to make the 
De-Broglie wavelength very short). Namely, in the case of 
a rough billiard there are eigenstates that have a lot of weight 
in the region \p\ < {2mEy/^. Consequently there are both 
semiclassical and non-semiclassical overlaps. Specifically, if 
we have non-semiclassical wavefunctions, and \En — E„i\ ^ 
0, then the collective contribution dominates, which give rise 
to the perturbative-like peak in the LDOS. 

Our findings apply to systems, such as the rough billiard, 
where there is an additional (large) non-universal energy scale 
i5i?Nu- This is defined as an energy scale which is not re- 



lated to the bandprofile, and hence does not emerge in the 
RMT modeling. Hence in general there is a distinct twilight 
regime h/ra < SE^,{(j)) < SE^^j, which is neither "perturba- 
tive" nor "semiclassical". [In our numerics ^=100 is so large 
that SEf,^,r^E.] 

Summary: We have analyzed the parametric evolution of 
the eigenstates of an Aharonov-Bohm cylindrical billiard, as 
the flux is changed. For the first time the full crossover 
from the perturbative to the non-perturbative regime is demon- 
strated. Random matrix theory suggests a simple crossover. 
Instead, we discover an intermediate twilight regime where 
perturbative and semiclassical features co-exist. This can be 
understood by adopting a phase space picture, and taking into 
account the inapplicability of the BeiTy conjecture regarding 
the semiclassical structure of the wavefunctions. 
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